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ABSTRACT 

Context. By studying young magnetically active late-type stars, i.e. analogues to the young Sun, we can draw conclusions on the 
evolution of the solar dynamo. 

Aims. We determine the topology of the surface magnetic field and study the relation between the magnetic field and cool photospheric 
spots in three young late-type stars. 

Methods. High-resolution spectropolarimetry of the targets was obtained with the HARPSpol instrument mounted at the ESO 3.6 m 
telescope. The signal-to-noise ratios of the Stokes /V measurements were boosted by combining the signal from a large number of 
spectroscopic absorption lines through the least squares deconvolution technique. Surface brightness and magnetic field maps were 
calculated using the Zeeman-Doppler imaging technique. 

Results. All three targets show clear signs of magnetic fields and cool spots. Only one of the targets, VI358 Ori, shows evidence 
of the dominance of non-axisymmetric modes. In two of the targets, the poloidal field is significantly stronger than the toroidal one, 
indicative of an o-^-type dynamo, in which convective turbulence effects dominate over the weak differential rotation. In two of the 
cases there is a slight anti-correlation between the cool spots and the strength of the radial magnetic field. However, even in these 
cases the correlation is much weaker than in the case of sunspots. 

Conclusions. The weak correlation between the measured radial magnetic field and cool spots may indicate a more complex magnetic 
field structure in the spots or spot groups involving mixed magnetic polarities. Comparison with a previously published magnetic field 
map shows that on one of the stars, HD 29615, the underlying magnetic field changed its polarity between 2009 and 2013. 
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1. Introduction 


Magnetic fields play a key role in the evolution of solar-type 
stars. They determine the angular momentum loss, shape stellar 
wind, produce high-energy electromagnetic and particle radia¬ 
tion, and influence the energy balance of planetary atmospheres. 

The general notion as postulated by iParked (1 19551) is that 
the solar magnetic field is generated by an crQ -dynamo: The 
toroidal field is generated by the winding up and simultaneous 
amplification of the poloidal field by differential rotation (Q ef¬ 
fect). The dynamo loop is closed by the collective inductive and 
diffusive effects arising from rotationally influenced convective 
turbulence, which generate small-scale poloidal fields from the 
underlying toroidal field (a effect), and the enhanced turbulent 
diffusion causing efficient reconnection to build up a large-scale 
poloidal field (j3 effect). In rapidly rotating late-type stars, the 
Q effect is suppressed and therefore other types of dynamos, 
name ly a^Q or a^, are probably present (see e.g. lOssendriiveil 
120031 and references therein). In this case, the magnetic field is 
almost solely sustained by the a effect generating both poloidal 
and toroidal fields. Typical ail solutions show dominance of the 
toroidal component, whereas the other types of solutions have 
more equal energies in both of the components of the surface 
field. 


Young late-type stars are generally rapidly rotating owing 
to the large angular momentum remaining from the contraction 
phase during early stellar evolution. Because of magnetic brak¬ 
ing the rotation slows down. This means that as a solar-type star 
evolves, its dynamo will move from the or a^Q -regime to 
an aO -dynamo. In order to properly understand the mechanism 
and evolution of the dynamos of solar-type stars it is essential to 
study stars of different ages. 

To investigate young solar-type stars, we initiated a pro¬ 
gramme called “Active Suns” for the observation a small sam¬ 
ple of young solar analogues with the HARPSpol instrument 
mounted at the ESO 3.6 m telescope at La Silla (Chile). The 
three stars in this study, AH Lep, HD 29615 and V1358 Ori, are 
located at the more rapidly rotating and thus more magnetically 
active end of our small sample. 

Our primary aims were to determine the topology of the 
surface magnetic field and to study the relation between mag¬ 
netic fields and cool spots. Eor this purpose we applied the 
Zeeman-Doppler imaging (ZDI) technique using Stokes /V 
spectropolarimetry with a code developed by O. Kochukhov 
(iKochukhov et al.ll2014l) . 


2. Observations 


* Based on observations made with the HARPSpol instrument on the 
ESO 3.6 m telescope at La Silla (Chile), under the program ID 091.D- 
0836. 


Stokes /V spectropolarimetry was collected with the HARPS¬ 
pol spectropolarimeter mounted at the ESO 3.6 m telescope at 
La Silla, Chile through the programme “Active Suns” (091.D- 
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0836). HARPS itself is a fibre-fed, cro ss-dispersed echelle spec¬ 
trograph described in greater detail by iMavor et al.l (l2003h . The 
pola risation unit HARPS pol is installed at the Cassegrain fo¬ 
cus dPiskunov et al.ll2()TTl) . The normal procedure for obtaining 
a Stokes V profile is to use four sub-exposures with a rotation of 
90° of the quarter-wave plate between each sub-exposure. This 
procedure allows the calculation of a “null spectrum”, which 
should only contain noise and is thus a way to check that there 
is no spurious signal contaminating the Stokes V profile. 

The observations were red uced using the REDUCE pack¬ 
age dPiskunov & Valenti|l2002h . The standard reduction proce¬ 
dure included bias subtraction, flat fielding by a normalised flat 
field spectrum, subtraction of scattered light, and an optimal ex¬ 
traction of the spectral orders. The blaze function of the echelle 
spectrometer was removed to first order by dividing the observed 
stellar spectra by a smoothed spectrum of the flat field lamp. 
The continuum normalisation was performed by fitting a second- 
order polynomial to the blaze-corrected spectra. A consistent 
continuum normalisation procedure was applied to the spectra 
extracted from all sub-exposures. The wavelength solution for 
each spectrum was accomplished by fitting a two-dimensional 
polynomial, for approximately 1000 extracted thorium lines ob¬ 
served from the internal lamp assembly. 

The spectral resolution of the reduced observations was R » 
105000. The rotation phase coverage (f^) for each star was es¬ 
timated assuming that each observatio n covers a phase range o f 
[^^rot-0.05, (^iot+0.05], as suggested bv iKochukhov et all(l2013h . 
The observations are summarised in Table[T]and listed in Table|2] 

Since the aim of the Stokes V observations was to map the 
surface magnetic field, the high demand for the signal-to-noise 
ratio (S /N) was far beyond what could be reached for individ¬ 
ual photospheric absorption lines within reasonable exposure 
times. In order to boost the S /N a large number of photospheric 
absorption lines in the wavelength range 3800 - 6900 A were 
combined with least squares deconvolution (LSD; iDonati et alJ 
Il997t iKochukhov et al.ll2010h . The line mask fo r the LSD pro¬ 
cedure was obtained from the VALD data base dPiskunov et al.l 
Il995t iKunka et alJll999ll . Regions dominated by strong lines 
were avoided, since the profile of strong lines deviates from the 
average line profile. More details on the LSD procedure can be 
found in the paper bv iKochukhov et ^ d20l4) . The final S/N of 
the LSD profiles was 5000 - 12000 in the Stokes V spectra. The 
LSD profiles were calculated on a radial velocity grid with a 1.6 
kms^' step size. The number of spectral lines used in the LSD 
mask for each star is listed in Table [T] The same masks were 
used for both the Stokes I and V. 

From the LSD Stokes V spectra we calculated the false alarm 
probabilities (FA P) for each observations using the red uced 
statistics (see e.g. lDonati et al.l[l992HRosen et alJl2013h and the 
radial velocity range [-vsin/-7.0 kms^*, vsin/-i-7.0 kms“^]. 
Each point of the Stokes V profile had its own error estimate cr, 
corresponding to the S /N. This cr, was used for the calculation 
of the FAP values. The same procedure was also followed for the 
null spectra, which naturally reached the same S /N level as the 
Stokes V. The minimum FAP values for the Stokes V profiles of 
the targets are listed in Table [T] and all FAP values for both the 
Stokes V and the null profiles are listed in Table|2] Using a FAP 
limit of 10“^ we concluded that all three targets showed definite 
detection of a magnetic field. We also evaluated the mean longi¬ 
tudinal magnetic field (B z) from the first moment of the Stokes 
V profile as described by IKochukhov et alJ (120101) for both the 
Stokes V observations and the null spectra using the same radial 
velocity range as for the FAP. These measurements are listed in 
Table |2] The {Bz) for each stellar observation as a function of 


Table 1. Summary of observations with HARPSpol at the ESO 3.6 m 
telescope. 


Star 

<5/V>“ 

wlsd* 

n/ 

// 

FAP ^ 

AHLep 

7400 

3493 

5 

0.49 

1.4- 10-'^ 

HD 29615 

9100 

3344 

9 

0.68 

< 10-'® 

VI358 Ori 

8800 

3085 

4 

0.40 

O 

1 

o 


Notes. Mean signal-to-noise ratio (S /N) of reduced Stokes V LSD 
line profiles, number of spectral lines in LSD (ulsd), number of 
observed phases (n^), estimated phase coverage (/^), minimum 
FAP for Stokes V signals. 

the rotation phase is plotted in Fig. [T] The errors on (Bz) were 
estimated by propagating the errors on the observations. The null 
spectra are plotted in Fig. |2 

3. Zeeman-Doppler imaging 

The Stokes V profiles were modelled assuming the weak field 
approximation, i.e. that the local Stokes V profile is proportional 
to the effective Lande g-factor. The Land e factor was determined 
as suggested bv iKochukhov et aTI(l2010l) and the values are listed 
in Table |3] 

The Stokes I profiles were modelled using a temperature in¬ 
dependent analytical Voigt function. For neutral photospheric 
metal lines, the absorption is usually larger in a cool spot than 
in the unspotted photosphere for solar-type stars. This effect is, 
however, of secondary importance because the reduced contin¬ 
uum level in the spot will dominate. Thus, for all photospheric 
absorption lines, the spot will be observable as an “emission” 
bump moving across the line as the star rotates. The amplitude 
of the bump will, of course, depend on the changing equivalent 
width of the absorption line. Therefore, the approximation of 
constant line profile should only lead to a slight error in the re¬ 
trieved brightness of the spots and the overall equivalent width 
of the line (as this usually becomes slightly stronger when large 
cool spots reduce the observed mean Tefr). 

Similar approximations to those used in this study are com¬ 
monly used in ZD I for the Stokes IV when applying the LSD 
technique (see e.g. lPetit et al.ll2004l) . An alternative would be to 
treat the LSD profiles as a single line with average parameters 
or derive the local LSD profile from the full spectrum synthesis. 
The former still involves approximations and so the advantages 
are limited. The latter is feasible for early-type stars (see e.g. 
IKochukhov et al.ir2014l) . but much more demanding for late-type 
stars when imaging both the surface magnetic field and temper¬ 
ature (or brightness). Furthermore, none of the stars in this study 
showed any strong bumps in the Stokes /, reducing the error of 
the approximation. 

We used the inversLSD code developed by O. Kochukhov. 
The stellar surface was divided into 1876 spatial elements. 
The magnetic inversion was bas ed on a spherica l harm onics 
expansion, a s desc ribed by e.g. IKochukhov et al.l (1201 4l) and 
iRosen etdl (l2015l) . One of the advantages with this approach 
is that one can easily determine the amount of energy contained 
in the toroidal and poloidal components of the magnetic field, as 
well as the different T’-modes. 

To set the highest allowed f-mode, different values of ^max 
from 3 to 10 were tested. In the solutions with > 5, smaller 
scale structures appeared even though the fits to the Stokes V ob¬ 
servations were not found to be significantly improved. This ten¬ 
dency was independent of the v sin /-value, which regulates the 
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Fig. 1. Mean longitudinal magnetic field as a function of rotational phase. From left to right: AH Lep, HD 29615 and V1358 Ori. 


Table 2. Observations with HARPSpol at the ESO 3.6 m telescope. 


Star 

Date 

f 

0rot 

^exp [s] 

SIN 

FAP [V] 

FAP [null y] 

(Bz) [G] 

Null (Bz) [G] 

AH Lep 

2013 Sep 9 

45.82 

0.918 

4x700 

8613 

1.4- 10-‘^ 

0.97 

-13.2 + 8.0 

-6.4 + 8.0 


2013 Sep 10 

46.81 

0.675 

4x700 

8822 

6.2 ■ 10-® 

0.72 

38.9 + 7.9 

-7.0+ 7.9 


2013 Sep 11 

47.81 

0.438 

4x700 

5265 

1.9-10-3 

0.94 

-20.6+12.7 

0.6+ 12.8 


2013 Sep 12 

48.81 

0.201 

4x700 

9744 

1.4- 10-3 

0.90 

18.9 + 7.2 

-8.6+ 7.2 


2013 Sep 15 

51.86 

0.523 

4x800 

4539 

0.44 

0.83 

12.8 + 14.7 

13.3+14.7 

HD 29615 

2013 Sep 9 

45.78 

0.837 

4x900 

10906 

< 10-^^ 

1.00 

-7.0+ 4.4 

1.5+4.4 


2013 Sep 9 

45.90 

0.887 

4x900 

9506 

< 10-'® 

0.81 

-11.1+5.1 

-3.2+ 5.1 


2013 Sep 10 

46.77 

0.264 

4x1000 

8524 

< 10-'® 

0.86 

-33.7 + 5.7 

1.1 +5.6 


2013 Sep 10 

46.89 

0.316 

4x1000 

10523 

< 10-'® 

0.91 

-25.4 + 4.6 

1.5+4.6 


2013 Sep 11 

47.77 

0.695 

4x1000 

7051 

0.32 

0.35 

-9.9+ 6.8 

-4.3 +6.8 


2013 Sep 11 

47.89 

0.746 

4x1000 

6025 

5.0 ■ 10-3 

0.26 

-0.0 + 7.9 

9.8+ 7.8 


2013 Sep 12 

48.77 

0.126 

4x1000 

12081 

< 10-'® 

0.84 

-38.3+4.0 

2.1 +4.0 


2013 Sep 12 

48.89 

0.177 

4x1000 

11542 

< 10-'® 

0.90 

-37.2+4.2 

4.1 +4.2 


2013 Sep 15 

51.81 

0.435 

4x1100 

5671 

< 10-'® 

0.94 

-13.4 + 8.4 

-7.6+ 8.4 

VI358 Ori 

2013 Sep 9 

45.86 

0.548 

4x675 

8994 

0.11 

0.72 

-14.7+13.4 

-5.4+13.4 


2013 Sep 10 

46.85 

0.403 

4x675 

9843 

1.7 ■ 10-'° 

0.98 

6.1 + 12.3 

-1.4+12.3 


2013 Sep 11 

47.85 

0.265 

4x675 

6390 

0.40 

1.00 

-7.1 + 18.5 

21.1 + 18.5 


2013 Sep 12 

48.85 

0.127 

4x675 

10050 

1.5 ■ 10-3 

0.91 

10.3 + 12.0 

9.7+ 11.9 


Notes. The time of each observation is given as HJD - 2456500 


expected surface resolution for temperature Doppler imaging. 
This implies that Stokes V alone is not sen sitive to sma l l-scale 
features of higher order f-modes. In fact, iRosen et"^ (1201 5h 
have shown that recovering such details would also require us¬ 
ing linear Stokes QU polarisation measurements. The maximum 
was therefore set to (max - 5 for all stars regardless of their rota¬ 
tional velocities. 

The magnetic field inversion was regularised with a har¬ 
monic penalty function, while we u sed Tikhonov regu larisation 
for the brightness inversion (see e.g. lRosen et al.llToiSh . In addi¬ 
tion to the Tikhonov regularisation we also introduced an extra 
penalty function to restrain the surface brightness within values 
of [0.5,1.5], with 1.0 corresponding to the unspotted surface: 

( ~ ^min) if bi < /imin 

f^p{b} ~ 1 ^ ^min — bi < ^max 

f X/ ^(^max “ bj^ if bi > /imax- 

Here bmin and b^n^y, are the allowed minimum and maximum val¬ 
ues for the brightness bi, and c is a constant. This extra restriction 
was necessary because using the approximation of a constant lo¬ 
cal Stokes I profile made the brightness inversion insensitive to 
the mean brightness. 


The inverse solution was obtained using the Levenberg- 
Marquardt minimisation. Here the first and second derivatives 
of the extra penalty function Rp were calculated analytically. 
The magnetic energy contained in the axisymmetric and non- 
axisymmetric harmo nic components w as calculated using the 
same definition as in iFares et al.l (l2009ll : Modes with m < (’/2 
were defined as axisymmetric and those with m > (jl consid¬ 
ered non-axisymmetric. 

4. Observed stars and their adopted stellar 
parameters 

The three stars in this study are nearby young solar-type and 
they are effectively single stars, in the sense that they lack any 
companions close enough for detectable interaction. AH Lep 
and VI358 Ori are believed to belong to the local association , 
implying ages of 20 - 150 Myr (see e.g. iMontes et al.l 120011) . 
while HD 29615 has been proposed as a member of the Tu- 
can a/Horologium associatio n indicating an age of ~30 Myr (see 
e.g. IZuckerman et'^l201 Ih . They are, however, slightly hotter 
than the young Sun. For one of the stars, VI358 Ori, there are 
no previously published ZDI studies. P revious ZDI maps o f AH 
Lep and HD 29615 were presented by ICarter et al.l (1201 5h and 
IWaite etal.1 (1201 .511 . 
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Fig. 2. Stokes V null profiles. The error bars are drawn as vertical lines 
through the data points and the 0-level is shown with a red line. From 
left to right: AH Lep, HD 29615, and V1358 Ori. 


Table 3. Summary of the adopted parameters. 


Star 

Pmt 

[d] 

V sin i 
[kms-'l 

R 

[Rs] 

i 

n 

Lande g-factor 

AH Lep 

1.31 

26.3 

1.21 

36 

1.214 

HD 29615 

2.32 

18.5 

0.96 

62 

1.214 

V1358 Ori 

1.16 

41.3 

1.1 

59 

1.207 


reported a photometric amplitude of Ay = 0.08. We estimated 
the rotation velocity to be vsin/ » 18.5 kms“* giving an incli¬ 
nation angle of i » 62°. Our observations were phased with the 
ephemeris 


HJD^=o = 2449730.0 H- 2.32 x E. 

Recently IWaite et ^ (1201 5h published a ZDI study of HD 
29615, where they also estimated the differential rotation of the 
star. Expressed as the differential rotation coefficient k - AQ/Q, 
they searched for the value of k that gave the best fit to the 
observations and reported k x 0.18 from the magnetic field 
structures and k 0.026 from the temperature spot struc¬ 
tures. In principle, a discrepancy between these values could be 
caused by different anchor depths. However, such a huge dif¬ 
ference seems unlikely. Further more, theoretical considerations 
(see e.g. [Krause & Radlei il98C ) and numerical sim ulations (e.g. 
iKitchatinov & Olemskov 1201 1 : ICole et al]l2014ll indicate that 
the magnetic fields of rapidly rotating solar-type stars are only 
weakly influenced by differential rotation, and they are therefore 
expected to have very small values of k. Thus, and especially 
since we simultaneously retrieve both the magnetic field and sur¬ 
face brightness, we did not apply any surface differential rotation 
in this study. 


4.1. AH Lep 


AH Lep (HD 36869, TY C 5916-792-1) is classified as a G3V 
star (iMontes e t al.l 200 1]) with a T YCHO catalogue parallax of 
28.6 mas (Wichmann et al.]l2003l) . i.e. a distance of r a; 35pc. 
IZuckerman et al.l ~i 201 ih reported an effective temperature of 
Teff ~ 5800 K and a radius of R ^ 1.21Ro. ICutispoto et al.l 
(1 19991) calculated a photometric rotation period of Pphot ~ 
lt'31. iMessina et al.l (l2003h reported a maximum photomet¬ 
ric amplitude of 0.09 in the V magnitude and an X-ray flux 
of lg(Lx/Lboi) = -3.480. Our best fit to the Stokes I pro¬ 
file was achieved with a rotation velocity of vsin/ = 26.3 
kms~*, i.e. close to t he va lue of 27.3 kms“* measured by 
iLopez-Santiago et al.l (l2010l) . Combining the previously mea¬ 
sured radius, rotation period and our rotation velocity estimate 
we retrieved an inclination of the rotation axis of / » 36°. To 
calculate the rotation phases (<prot in Tab.|2]i of the spectropolari- 
metric observations, we used the ephemeris 


HJD^=o = 2449730.0 -H 1.31 x £. 


4.2. HD 29615 


HD 29615 (HIP 2163 2, TYC 6467-702- 1) has a Simbad spectral 


classification of G3 V (Torres et 

alJl2006h. Its Hipparcos parallax 

is 17.8 mas (Ivan Leeuwenll2007 

) giving a distance of r » 56pc. 

According to McDonald et al.l ( 

2012^ its effective temperature 

is T’.-ff X, 5866X. Vidotto et al.| (|2014l) reported a radius of 


R ^ 0.96Rq and a rotation period of 2‘/32. IMessina et aP (l2010h 


4.3. VI358 Ori 


V1358 Ori (HD 43989, HIP 30030, TYC 4- 7 88-12 7 2-1) was clas¬ 
sified as a F9V-type star by iMontes et al.l (1200 ll) . Ivan LeeuwenI 
(l2007l) reported a Hipparc os parallax of 20 . 31 ma s, placing it 
at a distance of r ss 49pc. iMcDonald et al.l (1201 2l) reported an 


effective temperature of Tpaf x 6032 K, while the estimate of 

Vican & Schneidei 

(2014 

) was Teff ~ 6100 K. Furthermore, 

) reported a radius of » L05/?o, 

Vican & Schneidei 

(2014 


while IZuckerman et al.l (1201 ll) presented the value R » 1.08/?©. 
We adopted the value 1.1 Rq. From photometric observations 
ICutispoto et al.l (l2003h calculated a rotation period of Rphot ~ 


1^16 and estimated the photometric amplitude to be Ay = 0.08. 
Furthermore, they reported an X-ray flux of lg(Lx/Tboi) = 
-3.629. We estimated the rotation velocity to be vsin/ 41.3 
km s ', giving an inclination of / a; 59°. To phase our observa¬ 
tions, we used the ephemeris 


HJD^=o = 2449681.5 H-1.16 X £. 


5. Results 

The resulting ZDI-maps are presented in Figs. |3|5] and their 
corresponding Stokes I and V profiles in Fig. |6] The strongest 
magnetic field is detected on HD 29615 and the weakest field 
on V1358 Ori (Table |4|i. HD 29615 shows an almost axisym- 
metric structure with both the magnetic field and cool spots 
concentrated at high latitudes around the pole. In VI358 Ori 
the axisymmetric component of the field is stronger than the 
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Fig. 3. ZDI maps of AH Lep in equirectangular projection. The dashed vertical lines mark the observed rotation phases and the pointed horizontal 
line the limit of visibility due to the stellar inclination. 
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Fig. 4. Same as Fig.[3]for HD 29615 
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Fig. 5. Same as Fig.[3]for V1358 Ori 


non-axisymmetric. AH Lep, on the other hand, has clearly non- 
axisymmetric geometries of both the magnetic fields and cool 
spots. It should be noted, however, that the ZDI maps of AH Lep 
and VI358 Ori are based on poor phase coverages. 

The poloidal component dominates the surface helds of AH 
Lep and HD 29615, while the toroidal component is stronger for 
VI358 Ori. In terms of complexity, the { - 1 mode dominates 
the helds of HD 29615 and V1358 Ori, while ^ = 2 is strongest 
for AH Lep (containing 31.2% of the magnetic energy). The in¬ 
formation on the magnetic held of each star is summarised in Ta¬ 
ble |4] This table lists the maximum detected surface held (Bmax) 
and the percentage of magnetic energy contained in the different 
components of the held (toroidal/poloidal, axisymmetric/non- 
axisymmetric and f = 1 mode). 

We also calculated the linear correlation between the mag¬ 
netic held maps and surface brightness maps. The correlation be¬ 
tween the radial magnetic held (|Brl) and surface brightness {b) is 
especially interesting. In the solar case cool spots are, of course, 
concentrated in regions were strong magnetic helds penetrate the 
surface. This implies a strong negative correlation r(\Br\,b). We 
could only hnd a weak correlation for the three stars, and in the 
case of V1358 Ori the correlation coefhcient is positive, i.e. of 
an opposite sign than expected (Table|4l). Again, poor phase cov¬ 
erage may play a role. 


6. Conclusions 

We see considerable differences in the topology of the surface 
magnetic helds of the three stars in this study. Since the stars 
themselves are evolutionarily and structurally close, it is likely 
that these differences rehect temporal changes rather than differ¬ 
ent dynamo mechanisms. 

We only have three snapshots, which means that far reach¬ 
ing conclusions cannot be drawn from our study. However, the 
poloidal and toroidal components are of same order of magni¬ 
tude in all three target s, which is a strong indication of a dynamo 
of or type (cf. lKapvla et al.ll2013[ where dynamos of dif¬ 
ferent types are compared in direct numerical simulations). On 
the other hand, the clear dominance of the axisymmetric mode 
in two of the three targets L somewhat unexpected from mean- 
held dynamo theory (e.s. iKrause & Radleilll986h . Again, this 
could be due to the non-stationarity over time, i.e. caused by the 
targets being at different phases of the oscillatory dynamo cycle. 
This conclusions is not at odds with the proposed type of the 
dynamo, as e ven for dynamo s, oscillatory solutions can be 
expected (see lKanvla et al.ll2013L and references therein). 

The previous ZDI study of HD 29615 was b ased on ob¬ 
serva tions collected in November -December 2009 (IWaite et al.l 
|2015|) . These ZDI maps look quite similar to our maps. In both 
cases the magnetic held is dominated by a strong radial held 
centred near the rotational pole that covers the whole polar re¬ 
gion. Interestingly, the polarity in the 2009 map is opposite to 
that in our map. Thus, the star has undergone a polarity reversal 
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Fig. 6. The Stokes IV observations (plus signs/crosses) and profiles (solid line in red) calculated from the ZDI solutions. The error bars for the 
Stokes V observations are drawn as vertical lines through the data points. From left to right: AH Lep, HD 29615, and V1358 Ori. 


Table 4. Summary of magnetic field detection. 


Star 

^max 

[G] 

Poloidal 

[%] 

Toroidal 

[%] 

Axisymm. 

[%] 

Non-axisymm. 

[%] 

1 

[%] 

r(\BA,b) 

AH Lep 

133 

58.9 

41.1 

38.9 

61.1 

27.0 

-0.19 

HD 29615 

245 

62.8 

37.2 

87.6 

12.4 

77.2 

-0.26 

V1358 0ri 

111 

31.3 

68.7 

72.9 

27.1 

70.9 

0.21 


between December 2009 and September 2013. The brightness 
maps similarly show a high latitude spot structure clearly offset 
from the strongest radial field. 

When studying the correlation between the magnetic field 
and temperature spots, one must remember that the detectabil¬ 
ity of the field is weighted by the surface brightness. If this 
were the reason for the dark spots not coinciding with strong 
magnetic fields, one would expect the detected magnetic field 
strength to be correlated with the surface brightness. Only for 
the case of VI358 Ori do we observe such a weak trend, and 
even then this is most li kely biased by the poor pha se coyer- 
age. On the other hand, I Rosen & Kochukho^ (1201 2l) showed 
that the ZDI method is capable of reproducing both the mag¬ 
netic field and cool spot structures in cases where they coincide. 
In the case of our maps of AH Lep and especially of VI358 
Ori, the poor phase coyerage makes the comparison of mag¬ 
netic and brightness maps uncertain. How eyer, while some ex¬ 
ceptions exist (see e.g. lCarroll et al.ll2012h . in most cases where 
surface magnetic field and temperature, brightness, or spot occu¬ 


pancy maps haye been retrieyed from the same spectropolarimet- 
ric data, there seems to be a discrepancy bet ween the concentra¬ 
tion of cool spots and magnetic fields (see e.g. lJeffers et al.l201 lb 
iKochukhoy et aljr2013l : I Waite et ^l2015h . 

A possible explanation for the lack of correlation between 
cool spots and a detected magnetic field could be that the large 
spots seen in Doppler images are actually spot groups with 
mixed polarity. In this case the surface resolution of the ZDI 
is not able to distinguish the detailed magnetic field configura¬ 
tion. The images would thus only display a larger scale field. If 
this were the case, one would expect the lifetimes of the spot 
structures to be yery short as many magnetic flux tubes of oppo¬ 
site polarities packed tightly together would reconnect and decay 
rapidly. This would most likely manifest itself through line pro¬ 
file yariations, and the time scale of changes would be helpful 
in proyiding indirect eyidence of the topology of the unresolyed 
structures. Unfortunately, the limited time span of our obserya- 
tions does not allow for such an analysis. The situation is further 
complicated because temperature spots can be generated with- 
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out m agnetic fields t hroug h vo rtex instability , as de monstrated 
by e.g. Kapyla et alj ( 201 ll) and iMantere et alJ (1201 Ih . 

iRosen et*^ ( 2015h compared ZDI results retrieved using just 
Stokes IV with the full set of Stokes IVQU components for the 
RS CVn star II Peg. They found few differences in the tempera¬ 
ture maps, but significant differences in the magnetic field solu¬ 
tions. The retrieval of the meridional field component especially 
benefited from using the linear Stokes QU components. Further¬ 
more, using just Stokes IV resulted in an underestimation of the 
total magnetic energy. Therefore, we can expect that the real sur¬ 
face magnetic fields of our targets are stronger and much more 
complex than indicated by this study. 

Despite these shortcomings, the results of this study can be 
used to retrieve information on possible stellar magnetic cy¬ 
cles. Here it is especially interesting to follow the evolution of 
the poloidal and toroidal components and the complexity of the 
field, and to detect polarity reversals of the large-scale field. Such 
reversals have been seen in many stars and with this study we can 
add HD 29615 to this group. 
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